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ABSTRACT 

Non-local thermodynamic equilibrium (NLTE) effects in diagnostically important solar Eei lines are important due to the strong 
sensitivity of Ee i to ionizing UV radiation, which may lead to a considerable under-population of the Ee i levels in the solar atmosphere 
and, therefore, to a sizeable weakening of Ee i lines. Such NLTE effects may be intensified or weakened by horizontal radiative transfer 
(RT) in a three-dimensionally (3-D) structured atmosphere. We analyze the influence of horizontal RT on commonly used Ee i lines 
in a snapshot of a 3-D radiation magneto-hydrodynamic (MHD) simulation of a plage region. NLTE- and horizontal RT effects occur 
with considerable strength (up to 50% in line depth or equivalent width) in the analyzed snapshot. As they may have either sign and 
both signs occur with approximately the same frequency and strength, the net effects are small when considering spatially averaged 
quantities. The situation in the plage atmosphere turns out to be rather complex. Horizontal transfer leads to line-weakening relative 
to 1-D NLTE transfer near the boundaries of kG magnetic elements. Around the centers of these elements, however, we And an often 
significant line-strengthening. This behavior is in contrast to that expected from previous 3-D RT computations in idealized flux-tube 
models, which display only a line weakening. The origin of this unexpected behavior lies in the fact that magnetic elements are 
surrounded by dense and relatively cool down-flowing gas, which forms the walls of the magnetic elements. The continuum in these 
dense walls is often formed in colder gas than in the central part of the magnetic elements. Consequently, the central parts of the 
magnetic element experience a sub-average UV-irradiation leading to the observed 3-D NLTE line strengthening. 
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1. Introduction 


The Sun’s atmosphere is strongly structured at all spatial scales 
that can be resolved. In particular, the magnetic field in the 
Sun’s photospheric layers is concentrated a t scales tha t have 
only begun to be resolved r ecently (e.g. iLagg et aO l2010l: 

I Martinez Gonzalez et all l2012h . In the face of limited spatial 
resolution, spectral lines computed in 3-D model atmospheres 
resulting from radiation MHD simulations may be used to probe 
properties of magnetic features at a spatial resolution exceeding 
that of the best observations available today. However, the 
task is a difficult one and therefore simplifying assumptions 
have to be made. In the past, an often used simplification was 
the application of local thermodynamic equilibrium (LTE) to 
the calculation of atomic level populations. By doing so, the 
influence of the radiation field on the level populations and, 
therefore, on the line formation, was completely neglected. An 
important step towards a higher degree of realism in the calcu¬ 
lation of population numbers was the introduction of non-LTE 
(NLTE) calculations in one-dimensional (1-D) atmospheric 
models, still negl ecting the influence of horizontal radiative 
transfer (RT) (e.g.. | Athay & Liteslll972l:ISteenbock & Holweged 


1984 iBovarchuk et all II985I: ISolanki & Steenboch 
Thevenin & IdiartI 119991: IShchukina & Trujillo Buenc 


The inclusion of proper 3-D NLTE RT reflects the next step in 
improving the realism of line profile calculations. 
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In iHolzreuter & SolanO (1201 2h (hereafter referred to as 
Paper I), we investigated the effects of 3-D NLTE RT in 
several idealized (i.e. thin-tube) models of flux sheets (ES) 
and flux tubes (ET), qualitative ly confirming the results of 
IStenholm & Stenflol (I1977L ri97^h who found that the line pro¬ 
files of iron lines in a sim plified ET may be strongly weak- 
ened by horizontal RT (cf. iBruls & von der Luh3 l2Q0ll) . The 
strong UV irradiation from the surrounding hot walls of the 
ET represented the main factor producing this l ine weakening 
by the so-called UV over-ionization (see, e.g., lAthav & Litesf 


Il972l:l^ttenlll988l: [Shchukina & Trujillo Buenoll200lL and ref- 

erences therein). However, the true efeect was found to be much 
smaller than originally calculated bv IStenholm & Stenflol (1 19771) 
because the UV opacity used in their calculation was too low as 
they neglected the opacity contributions from the uncounted iron 
absorption lines in the UV wavelength range. 


A natural next step towy ds increased realism was taken by 
IHolzreuter & SolanO (1201 3l) (hereafter called Paper II), where 
we used a snapshot of a 3-D radiation hydrod ynamic (HD) 
simul ation computed with the MURAM code (IVogler et al.l 
20051) as the atmospheric in put, thus extending the work of 
Kiselman & Nordlimdl (1 19951) . In such an atmosphere, the line 
weakening effect by 3-D NLTE RT due to hot UV irradiation, 
acting mainly in intergranular lanes, occurs alongside the oppo¬ 
site effect of line strengthening in regions with colder surround¬ 
ings, i.e. mostly at the centers of granules. Such a strengthening 
was found already in Paper I at a few positions, but with a negli- 
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gible influence on the spatially averaged line profile. The much 
stronger and more frequent occurrence in the HD atmosphere 
leads to very similar average profiles of the 1-D and 3-D NLTE 
calculations. The spatially resolved line profiles, however, may 
differ considerably between the three calculation methods (LTE, 
1-D NLTE, and, 3-D NLTE) with corresponding implications for 
the interpretation and inversion of high-resolution observations. 

Here, we continue our investigation on the effects of true 
3-D RT on diagnostically important iron lines by considering 
a snapshot of a magneto-hydrodynamic (MHD) simulation of 
a plage region calculated by the MURAM code. The introduc¬ 
tion of a (strong) magnetic held in the MHD simulation leads to 
the formation of evacuated magnetic elements. We show that th e 
heavily idealized models of [Stenholm & Stenflol (I1977L Il978l) : 
iBruls & von der Liiha (l200ll) and Paper I can only partially re¬ 
produce the far more complex situation of magnetic elements 
in an MHD atmosphere. Note that the 3-D formation of chro¬ 
mospheric spectral li nes in 3-D radiat i on MHD simulatio n has 
been investigate d by iLeenaarts et al.l (l2009l 120121 l2013l) and 
iLeenaarti (l201Ql) 


2. Model ingredients 

2.1. Radiative transfer calculations 

All calculations were done in the same way as in Paper II, ex¬ 
cept that the full Stokes vector of the Zeeman split lines was 
computed. Eor a detailed description, especially of the formal 
solution of the full Stokes vector problem, we refer to Paper I 
(Appendix I and II). Here, we present only a short summary. 

We used the RH code (IUitenbroekl200Ql) with adaptations as 
mentioned in Papers I and II, the most important being the inter¬ 
polation method of the formal solution along the short character¬ 
istics, where we used mo notonic parabolic Bezier integration as 
proposed by lAued (l2003l) together with an enhancement we in¬ 
troduced for the magnetic case (see the Appendix B of Paper I). 

The atomic model was the same as used in Paper II (23 
levels, 33 lines). It was tailored to optimally reproduce the 
NLTE effects on the population numbers of levels involved in 
the formation of the transitions investigated here, i.e. the Eei 
524.71 nm/525.02 nm and the 630.15 nm / 630.25 nm line pairs, 
and, at the same time, to be as small as possible to reduce the 
computational effort. The pro blem of strong UV over-ionization 
owing to missing UV opacity (iBruls et al.lll992h was — again — 
eliminated by artificially increasing the opacities in the relevant 
wavele ngth range (opacity fudging) as proposed by IBruls et ^ 
(Il992h . Refer to Paper I for a more complete discussion. 

The full Stokes vector was calculated in the same w ay as 
in Pap er I, i.e. by using the so-called field-free method of iReesI 
(1 19691) . Thus, the influence of the Zeeman effect on the line pro¬ 
files (splitting, polarization) was included only in the final deter¬ 
mination of the emerging intensity but not during the iteration 
of the level population numbers. In the following, however, we 
analyze only the Stokes I profiles since we are interested here 
mainly in uncovering the basic effects. The influence of neglect¬ 
ing 3-D NLTE on inversions of Stokes profiles will be considered 
in a following paper. 

2.2. Model atmospheres 

The input model for our calculations was taken from a realistic 
3-D radiation magneto-hydrodynamic simulation of a plage re¬ 
gion. We used a snapshot from slB = 400 G run of the MURAM 


code (IVogler et al.ll2005l) . i.e., a run with an initial unipolar (pos¬ 
itive polarity) vertical homogeneous magnetic held that was im¬ 
posed on a relaxed HD run. The original resolution of the cube 
was 288x288x 120 voxels spanning a geometric range of ap¬ 
proximately 6 Mm X 6 Mm X 1.6 Mm. 

The construction of the atmospheres (3-D and 1-D) was per¬ 
formed in the same way as in Paper II. However, only the bottom 
30 grid points were removed instead of 48 as in the HD case of 
Paper II) owing to the strong fields in the present atmosphere and 
the associated evacuation, which pushes the log r = 0 surface to 
deeper layers. The resulting 288x288x90 cube starts approx¬ 
imately 500 km below the average = 1 level (z = 0 km) 
and reaches a height of approximately 750 km above it. The 
minimum value at the bottom of the reduced cube is ap¬ 
proximately 17, while the average value at that depth amounts to 
above 1000. 

Maps of the most important physical parameters of the in¬ 
put atmosphere are given in Pig. [T] for the same two geometric 
heights as in Paper II, i.e. at z = 0 km and z = 240 km, the latter 
very roughly being the height of formation of the cores of the 
lines used in this work. Note that the granular-intergranular tem¬ 
perature inversion also takes place at these heights. Some parts 
— mostly in the lower half of the image — still show hot up¬ 
rising granules, whereas those in the upper part consist of hotter 
intergranular lanes. The panel in Pig. [T] showing the vertical ve¬ 
locity (v^) confirms this view, as the granules in the lower part 
still show considerable upward velocities (blue), whereas those 
in the upper part of the image barely host rising gas at all (grey). 
Note also, that the downward velocities are often higher at the 
boundaries of the magnetic elements, whereas the velocities in 
the centers of the magnetic elements are often reduced owing to 
the magnetic held, although at some locations very strong down¬ 
flows are found inside the magnetic features. 

Maps of the vertical component of the magnetic held, (top 
panels of Pig. [T]), reveal held concentrations along intergranu¬ 
lar lanes (in the form of strongly elongated magnetic elements 
similar to flux sheets) and at the intersections of intergranular 
lanes (similar to flux tubes). Eor simplicity in the following we 
refer to the elongated magnetic elements as flux sheets, or PS, 
and to the rounder ones as flux tubes, or ET. Only few areas are 
found with negative vertical fields and the held strength there is 
much lower than that in the magnetic elements. They are mostly 
located at the boundaries between intergranu lar and granular ar¬ 
eas where th e down-flows are strongest fsee lBiihler et al.|[^Q15l: 
lRiihleHl2ni4 . 


3. Results 

3.1. Continuum intensity 

The computed intensity at a continuum wavelength close to 
525 nm is presented in Pig. [21 Two types of regions of more in¬ 
tense radiation can be found: One in the hot granules and one in 
the cores of magnetic flux concentrations. Lower intensities are 
present in the intergranular lanes between magnetic PS and gran¬ 
ules. The results of our three calculation methods (3-D NLTE, 
1-D NLTE, LTE) show no differences in the continuum intensity 
as it is formed strictly in LTE in our computation. Note that we 
neglec t the possible departures from LTE found bv IShapiro et al.l 
(|201Q|) . The two black lines indicate cuts through an ET and an 
PS, respectively, which are investigated in detail in later sections. 
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Fig. 2. Continuum intensity at 525 nm. The grey scale has 
been clipped to 10% and 90%, respectively, of the maximum 
continuum intensity value. The two black lines indicate the se¬ 
lected cuts through two distinct magnetic elements, namely, an 
elongated feature similar to a flux sheet (FS) and a more compact 
one that we will call a flux tube (FT). 


3.2. Influence of NLTE and horizontal RT 

In this section, we compare the results of LTE, 1-D NLTE 
and 3-D NLTE radiative transfer on the spectral lines. We 
are looking for the effects of ho rizontal RT as described by 
IStenholm & Stenfld (Il977l Il978l) . who found a pronounced 
weakening of the lines calculated in 3-D NLTE in magnetic el¬ 
ements. Such a weakening may lead to errors when diagnosing 
atmospheric parameters such as the temperature. To simplify the 
discussion, we introduce the following notation for the relative 
differences of intensities (I) of the different calculation methods 
(3-D NLTE, 1-D NLTE and LTE): 





( 1 ) 


with each M/, and Mj referring to one out of the three cal¬ 
culation methods (3D, ID, LTE), and being the lowest in¬ 
tensity value of the line obtained with method M/ at a specific 
spatial position. In the following, we will call residual inten¬ 
sity of method Mf. Analogously, for the equivalent widths (EW) 
we define 




- EW^j 
E^j 


( 2 ) 


Note that > 0 and < 0 if the 3-D NLTE 

line is weaker than its 1-D NLTE counterpart, i.e., is a 

measure for the line weakening of the 3-D NLTE line when com¬ 
pared with the 1-D NLTE line strength. The larger the 
is, the stronger is the line weakening. Negative denote 

locations with line strengthening. 


3.2.1. Equivalent width 

Eigure [3] shows the spatial distribution of (left pan¬ 
els) and (right panels). Upper panels: 525.02 nm; lower 

panels: 630.15 nm. Yellow and red colors indicate areas where 
the 1-D NLTE (3-D NLTE) line is weaker than the correspond¬ 
ing LTE (1-D NLTE) line. The results for the 524.71 nm and 
the 630.25 nm lines are not shown, since they agree qualitatively 
with those of the other line in the respective line pair. The effects 
in the 630.15 nm line are slightly weaker than those found in the 
630.25 nm line. 

Looking at Eig. [3l we And large values in all panels. The 
influence of NLTE (left panels) is typically larger than that of 
horizontal RT (right panels) and often has the reversed sign, i.e., 
NLTE effects are partially cancelled by those of horizontal RT if 
one compares 3-D NLTE with LTE calculations. 

Where are these differences located spatially? Eirst consider 
the distribution of in the left panels of Eig. O In the 

granules and in the centers of the magnetic elements we And 
^ELTe-id ^ whereas < 0 in the periphery of 

the magnetic elements. The positive and negative deviations ap¬ 
proximately reach the same strength (up to 20%) for the Eei 
630.15 nm line. Eor the 525.02 nm line, is very large in 

the granules and the centers of the magnetic elements and close 
to zero elsewhere. 

This behavior results from the fact that the relatively weak 
525.02 nm line’s source function (LTE and NLTE) mainly fol¬ 
lows the local Planck function. At locations with a steep vertical 
temperature gradient (e.g. at a granule), UV irradiation from be¬ 
low ionizes iron and thus reduces the population of the lower 
level. The resulting lower opacity of the NLTE line shifts its for¬ 
mation height towards deeper and hotter regions in the atmo¬ 
sphere, leading to an increase of the source function and there¬ 
fore a weakening of the line in NLTE. In regions with smaller 
temperature gradient, the effect of NLTE disappears almost com¬ 
pletely. These arguments apply to the Eei 630.15 nm line as 
well. However, for the Ee 1 630.15 nm line an additional effect is 
found that acts with a reversed sign. The 630.15 nm line is much 
stronger than the 525.02 nm line. Therefore, it is formed higher 
and the source function falls clearly below the Planck function, 
making the NLTE line generally stronger. This balance between 
the two effects (over-ionization and low source function) leads 
to the pattern seen in the lower left panel of Eig.[3l with the line 
being stronger in LTE in the granules, but weaker in the inter¬ 
granular lanes. 

The right panels of Eig. [3] display the influence of horizontal 
RT. Both lines exhibit almost the same pattern, the line weak¬ 
ening in 3-D NLTE being slightly more pronounced for the Eei 
630.15 nm line. In the magnetic elements, with the exception 
of their central parts, the 3-D NLTE lines are weaker (red and 
yellow parts in the right panels of Eig. [3]), a s expected from the 
investigation of IStenholm & Stenflol (11^71) and Paper I, owing 
to the strong irradiation from the hot environment. 

However, the situation in the central parts of the magnetic el¬ 
ements is reversed (small blue and green patches within the red 
and yellow lanes in the right panels of Pig. [3]) : At these locations, 
we And the 3-D NLTE lines to be clearly stronger. This is a sur- 
prisin g result that contradicts the r esult from the idealized mod- 
els of istenholm & Stenflol (ll977l) : lBruls & von der Liihd (l200lh 
and Paper I, where the line weakening effect of horizontal trans¬ 
fer was strongest in the centers of the elements. Obviously, the 
idealized models differ significantly from the more realistic pic¬ 
ture provided by the MHD simulations. 
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In large parts of the granules is negative, a result al¬ 

ready found in Paper 11. The origin of this effect lies in the influ¬ 
ence of the cooler horizontal environment near the line forming 
region, i.e. the intergranular lanes. 

Figure |4] depicts histograms of the weakening/strengthening 
of the 525.02 nm and the 630.15 nm lines calculated in LTE and 
1-D NLTE in relation to the 3-D NLTE line profile. While the 
top panels present the distributions for the whole snapshot, the 
middle and lower panels shows histograms for selected subsets 
of the atmosphere. In the middle panels, the selection represents 
areas of strong up- (thin) and down-flows (thick), respectively. 
The selected points all have |v^| > 3 km/s at the Tc = I level. 
As expected from Fig.[3l we And > 0 in the down-flow 

areas and < 0 in the up-flows areas and the reversed sit¬ 
uation for the because the LTE - 1-D NLTE difference 

is stronger than that between 1-D and 3-D NLTE. The width of 
the distribution is often quite large, especially for the 

525.02 nm line. 

In the bottom panels we distinguish between the subset of 
pixels with > 1 kG (called ’’magnetic”) and those with 
\Bz\ <0.1 kG (called ”non-magnetic”). Interestingly, we obtain 
that the average line strengthening and weakening is almost the 
same for the magnetic and the non-magnetic areas. However, 
the width of the distribution is much larger for the magnetic 
areas, most magnetic pixels are associated with a line weaken¬ 
ing in 3-D, but the distribution has a long tail towards great line 
strengthening. Strong down-flows which mainly take place out¬ 
side magnetic elements harbor line weakening in 3-D NLTE as 
they have hot surroundings. Furthermore, as mentioned above, 
in the innermost parts of magnetic elements considerable line 
strengthening takes place. 

3.3. Line formation in a realistic flux tube 

In the previous section we found lines formed in 3-D NLTE to be 
stronger than their 1-D NLTE counterparts in the centers of mag¬ 
netic elements (i.e. < 0, in contradiction to previous 

results from idealized models. In order to better understand this 
surprising result and to compare it with the results of Paper I, we 
selected a sample cut through a magnetic element that is FT like. 
The cut runs diagonally and is centered onx^^ L8,yj^3.3 Mm 
as indicated in Fig.[2l 

3.3.1. Line profiles 

Profiles of Fei 525.02 nm and 630.15 nm are plotted in Fig. [5] 
at five spatial locations along the cut through the FT used in the 
previous section (see Fig. 0. The profiles between these posi¬ 
tions typically show a continuous, slow change towards the next 
plotted profile. At position x = 2.25 Mm, in the middle of a gran¬ 
ule, the profiles are deep, with a strong continuum. Small differ¬ 
ences on the order of 1 % to 3% (LTE) are found between the line 
depths of the three profiles at this spatial location, the 1-D NLTE 
line being slightly weaker than the 3-D NLTE line, as expected 
from Paper IT At position x = 1.85 Mm near the boundary of the 
ET, the line intensities as well as the continuum are much weaker 
and, owing to the Zeeman effect, the line is broader. At this po¬ 
sition, we And > 0, i.e. line weakening, in agreement 

with idealized ET. The Ee i 525.02 nm line, however, is so weak 
and formed so deep that only minor absolute differences remain. 
Note the different y-axis ranges. At x = 1.77 Mm, in the inner 
ET, the lines are even broader and the Ee 1 525.02 nm line is com¬ 
pletely split by the Zeeman effect. Here we And ^ 0. 
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Fig. 6. Two extreme examples of the influence of NLTE and 
horizontal RT on line profiles: left panels: Eei 525.02 nm line; 
right panels: Eei 630.15 nm line; upper panels: x 2.15 Mm, 
y 0.4 Mm; lower panels: x 3.61 Mm, y = 1.50 Mm. Line 
styles as in Eig.[5] 


In the core of the ET, at x = 1.70 Mm, is negative 

for the 630.15 nm line, and the 525.02 nm line is so weak that 
no apparent difference may be seen. At x = 1.64 Mm we And 
~ 0 on the other side of the ET. At spatial locations be¬ 
tween X = 1.64 Mm and the outer boundary of the ET 
is again positive, representing the ’’normally expected” case of 
line weakening. 

The profiles along our cut represent standard cases. The dif¬ 
ferences between the individual profiles are mostly small. There 
are also locations within the analyzed snapshot with much larger 
differences. Eig. [6] presents two such cases where \6E^^~^^\ 
reaches up to several 10%. Although the differences in the pro¬ 
files are much larger, the two locations still represent analo¬ 
gous situations: The profiles in the upper panels originate from 
the center of a ET, whereas those shown in the lower panels 
emerge from a ET boundary. Note, that the LTE profiles are 
much stronger in both cases, i.e. the inclusion of horizontal RT 
may either increase or decrease the differences to LTE values. 

3.3.2. Atmospheric conditions in the sample flux tube-like 
structure 

Eigure|7]represents the temperature profile (upper panel) and the 
vertical magnetic held strength (lower panel) along the cut. The 
dotted lines in white indicates the continuum (jc = 1) level. At 
most horizontal locations (granules, center of ET/ES) the con¬ 
tinuum is formed in the temperature range above 6300 K, i.e. 
the region colored in red. However, it is formed in a cooler layer 
(green layer with 6000 K) close to the borders of the ET where 
the vertical gradient of the temperature is smaller (compare the 
thickness of the green and red layers at the different positions, 
e.g., at X = 1.8 Mm with x = 2.2 Mm). This corresponds well 
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Distribution of the vertical field 



Fig. 7. Temperature distribution (upper panel) and vertical mag¬ 
netic field Bz (lower panel) along the cut through the FT indi¬ 
cated in Fig. [2] (diagonal line). The colored dots indicate several 
T = 1 levels (white: continuum near 525 nm; yellow: 525.02 nm 
line core; brown: 630.15 nm line core). Note that the values of 
the x-axis correspond to those of the x-axis of the simulation and 
not to the length along the diagonal cut. I.e., the true distance be¬ 
tween two points along the (diagonal) cut is longer by a factor 

V2. 



1.4 1.6 1.8 2.0 2.2 

X [Mm] 


Fig. 8. Physical quantities at a fixed geometrical height (z = 
-42 km) along the same diagonal cut as in Fig. [71 This height 
corresponds roughly to the line core formation height of the 
630.15 nm line in the center of the FT. The axis description 
of the solid lines is given on the left, that of the dashed lines 
on the right side. Top panel: Temperature (solid) and density 
(dashed). Middle panel: Magnetic field (solid) and vertical ve¬ 
locity (dashed, positive values correspond to upfiows). The tick 
marks at the horizontal zero line correspond to the individual 
voxels in the model atmosphere. Bottom panel: NLTE intensity 
ratio (/ 3 _p) - ^1 -d)/^ 3-D (solid) representing 3-D NLTE line 
weakening (positive values) or strengthening (negative values) 
and continuum intensity (dashed). Vertical lines (dotted) repre¬ 
sent positions (numbered on top of the uppermost panel) along 
the cut that are discussed in the main text. As in Fig. [71 the val¬ 
ues of the x-axis correspond to x-axis of the simulation and not 
to the length along the diagonal cut. 


with the brightness of the continuum in Fig.[2l By comparing the 
two panels of Fig. [ 7 ] we clearly see that the r = 1 level starts to 
curve down well outside the magnetic feature. Only a small part 
of this change in r = 1 outside the FT is due to the expansion 
of the FT, which places low density gas above these locations 
(see Paper I). It is mainly due to the lower temperature in the 
downflow lanes. 

The yellow and brown lines in Fig. [ 7 ] denote the r = 1 level 
of the 525.02 nm and the 630.15 nm line cores, respectively. Due 
to its strong temperature sensitivity, the 525.02 nm line is formed 
almost at the same height as the continuum over large parts of 


the FT thus minimizing the possible influence of NLTE effects 
in general. The height of its formation in the granules amounts 
to about 200 km above = 1, i.e., the Wilson depression in 
the 525.02 nm line core is far larger than that in the continuum, 
whereas it is only about 100 km larger than that of the continuum 
in the 630.15 nm line core. 

Fig. [8] presents some quantities along the cut through the FT 
at a fixed geometrical height (z = -42 km). The 630.15 nm line 
core is formed at this approximate height in the FT. For bet¬ 
ter visibility and comparison of the different panels, some of 
the horizontal x-positions that are mentioned in the following 
are marked with a thin dotted vertical line (numbered on top of 
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the uppermost panel). Several interesting facts can be seen from 

Fig.® 

a top panel: The drop in temperature (thick line) associated 
with the FT starts at x = 1.40 Mm (pos. 1) and ends at 
X = 1.92 Mm (pos. 5). This drop is seen at a fixed geomet¬ 
ric height. There are two boundary regions with enhanced 
density (dashed line) in the cool part of the FT of the atmo¬ 
sphere. 

b top and middle panels: The maximum density is reached at 
the positions 2 and 4 just outside the boundary of the FT 
where the magnetic field (thick line) has dropped to nearly 
zero. 

c top and middle panel: According to the magnetic field 
strength the FT diameter at the chosen height amounts to 
approximately 470 km (taking into account the V2 factor 
because the cut is in diagonal direction). The region with 
low temperature is considerably larger, i.e., approximately 
650 km. 

d middle panel: The drop in temperature occurs over a small 
range in x at the boundary between granular up- and down¬ 
flows, i.e, at the edges of the granules. 

e bottom panel: The spectral line formed in 3-D NLTE is 
weakened relative to that in 1-D NLTE (thick line, posi¬ 
tive values) in most locations along the cut, except in the 
strongest field part of the FT and in the granule to the right 
(x > 2.1 Mm). The strengthening in the granule is expected 
from the results of Paper II). The granule on the left is small 
and has a complex structure (see continuum image in Fig.O, 
so that the effect is not so clean. 

f bottom panel: The 3-D NLTE line is strengthened in the core 
of the FT (pos. 3). The temperature is also slightly higher at 
that location. 

g bottom panel: The continuum intensity (thin line) roughly 
corresponds inversely to the line weakening/strengthening. 
The locations of lowest continuum intensity coincide with 
peaks in the gas density. 

According to Fig. 1 in Paper I the continuum in the ideal¬ 
ized FT is formed at a depth at which the surrounding walls are 
at a temperature of 10000 K. In contrast, the walls in the MHD 
atmosphere, i.e. the region where the density increases and the 
field strength decreases, harbor down-flows which are not sig¬ 
nificantly hotter than the interior of the magnetic element at the 
same geometric height. This cool gas immediately surrounding 
the FT has an even higher density than that in the granules at 
the same geometrical height. The higher temperature reached in 
granules is located further away from the magnetic feature, be¬ 
hind the humps of high density. 

The 3-D NLTE line strengthening in the core of the ET is 
a consequence of the walls being not as hot as in the idealized 
ET case. Contributions from the very hot granular gas (with T > 
7000 K) are absorbed before they reach the inner part of the 
ET. Eurthermore, the immediate surrounding of the ET core is 
cooler on average than the core itself. This is illustrated in Eig.® 
where the temperature curve at the ET core (bold line) and the 
run of the average temperature of the surroundings at different 
distances and directions (thin lines, see figure caption) is given. 
The temperature in the ET core at the height of formation of the 
630.15 nm line core (vertical line) is several 100 K higher than 
that in its immediate surrounding. The run of the temperature 
much swallower than in the atmosphere outside the ET. These 
two facts lead to the situation that an atom at the center of the ET 
senses less UV irradiation than an atom closer to the boundaries. 



z [km] 

Fig. 9. Temperature curves as a function of height in the core 
of the chosen ET where the 3-D NLTE line strengthening is 
strongest (thick solid line) and at positions in the surroundings of 
the ET. Other solid lines: temperature curves averaged over the 
four points in positive and negative x and y direction at a given 
distance from the ET center. The thickness of the line increases 
with distance from the ET core. Erom thinnest to thicker lines 
(solid): 84 km, 146 km, and, 209 km. Dashed lines: the same as 
for the solid ones but for the four points in diagonal direction. 
The two lines lie 118 km (thinner line) and 177 km (thicker line) 
away. The vertical line indicates the approximate height of for¬ 
mation of the 630.15 nm line core in the core of the ET, where 
the continuum is formed slightly below z = -200 km. 


3.4. Conditions in a weak flux sheet 


The MHD snapshot used in this work, shows two typical forms 
of ES. The two forms can be distinguished in the continuum im¬ 
age of Eig. [2l Many intergranular lanes are dark and have a cen¬ 
tral brightening. They show a similar behavior like that of the 
ET discussed in the previous section. The second form of inter¬ 
granular lanes, which will be discussed here, shows no central 
brightening in the continuum. Eigure[T0| presents several atmo¬ 
spheric quantities along a cut through such a ES, which, in this 
case, is tilted towards a neighboring granule, a situation quite 
common in the analyzed simulation snapshot. The main reason 
for the difference in behavior is the lower field strength in this 
ES (see Tritschler & Uitenbroekl 120061: iRiethmliller et al.|[20l4 
iRiethmuller & Solankill2015li rThe temperature in the ES (left 
panel) is quite similar to that in the ET. However, the magnitude 
of the temperature fluctuations is smaller, as expected from the 
weaker magnetic field (see middle panel). In the left panel the 
continuum and line formation height is indicated, as in Eig. [T] 
The continuum is formed in two different layers. In the granules, 
the formation height is in the temperature range above 6300 K 
(red parts), whereas it is within layers of approximately 5600 K 
(green layer) in the down-flowing parts of the atmosphere closer 
to and in the ES. 


The 525.02 nm line displays again a much larger Wilson de¬ 
pression than the continuum. The Wilson depression in the core 
of Ee 1 630.15 nm is of intermediate magnitude. Owing to the tilt 
of the ES, the location of deepest origin follows the location of 
the strongest field at the height of formation, shifting from the 
left to the right in the left panel of Eig. [TT] when going from the 
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very close to zero near the middle of the FS. If one compares 
with Fig.[7]and[3]one finds that the line weakening is strongest in 
the intergranular lanes but outside magnetic features. It is small¬ 
est in the magnetic features, irrespective of whether the FS core 
is bright or dark in the continuum (Fig. O. There are only a 
few exceptions to that rule. E.g. the lane near (x = 2.6 Mm, 
y = 2.2 Mm), for which the line weakening is most pronounced 
in the middle of the FS. However, this FS has B <= 500 G at all 
heights in the photosphere. Such small B values lead to a minute 
Wilson depression so that B does not influence the horizontal RT 
in a noticeable manner. 

In summary, we find that flux sheets display line weakening 
du e to horizontal irradiation i n magnetic features as described 
bv IStenholm & Stenfl^ (Il977l) . although in a mild form, due to 
the presence of dense cool material on both sides of the FS. The 
inverse effect found in the centers of the FT, is less pronounced 
or not present owing to the lower field strength. 



y [Mm] 

Fig. 11. The same quantities as in Fig. [8] along the cut through 
the FS at geometrical height z = 0 km. 


630.15 nm towards the more weakened 525.02 nm line and the 
continuum. 

From the middle panel of Fig.[T0l we see that B^ is of the or¬ 
der of several tiG in the FS. The maximum value amounts to ap¬ 
proximately 1 kG. Therefore, the corresponding Wilson depres¬ 
sion is much smaller than that in the FT discussed in Sect. 13.31 
The FS is surrounded by two regions of reversed B^ of weaker 
field strength (150 G, dark blue^ The strong magnetic field in 
the FS does not extend as deep as that of the FT. The strong tur¬ 
bulence of the granular atmosphere deforms the relatively weak 
field and breaks up the sheet in deeper layers. This is typical of 
many FS-like structures and the one shown here is representative 
in this respect. 

Figure [TT] presents the same quantities at height z = 0 km 
as Fig. [HI but now along the cut through the FS. The situation is 
very similar: There are again two cool, density enhanced walls 
on both sides of the FS. The density decreases as the magnetic 
field increases, but interestingly remains above the value in the 
granules on either side. The width of the magnetic FS is much 
smaller than that of the region of cool gas, i.e. the walls are, 
at a given height, mainly at the same temperature as the inte¬ 
rior of the FS. The strongest down-flow is found at the boundary 
of the FS where the magnetic field direction is reversed (proba¬ 
bly due to the strong down-flow pulling the weak magnetic field 
downwards). Such reversed polarity fields s urrounding magneti c 
elements were d etected obseryationally by IZaver et M d 19891) : 
iNaravanI (|201 ll) : [Bfihleil (l20l4 : lBuhler et all (1201^ 

From the bottom panel of Fig.[TT]we see also that line weak¬ 
ening > 0) takes place over the full width of the FS and 

in its immediate surroundings. However, dips to values 


4. Discussion 

In their pioneering investiga tions of 3-D radiati v e transfer in 
heavily idealized flux tubes, IStenholm & Stenfl^ (1 19771 Il978l) 
found considerable net weakening of Fei lines owing to the 
strong horizontal UV-irradia tion from th e hot walls of the 
strongly evacuated flux t ubes (l^ruitll 1 976l) . Thi s result was cor¬ 
roborated by the work of lBruls & von der Liitiei (1200 ll) in ideal¬ 
ized 2-D and, in Paper I, in idealized 3-D models of FS and FT, 
although the effect was clearly weaker in these newer studies be¬ 
cause of the more realistic UV opacity used. In Paper II, using a 
snapshot of a HD simulation as the input atmosphere, we found 
line strengthening, owing to below-average UV irradiation from 
colder environments, to be as frequent as line weakening due to 
excess UV irradiation from hot surroundings. To investigate the 
situation in a plage region, we here applied our full 3-D NLTE 
RT calculation to a snaps hot of an MHD sim ulation carried out 
with the MURAM code (IVogler et al .1120051) with an initial ho¬ 
mogeneous vertical field of 400 Gauss. Naively one would ex¬ 
pect the line weakening to dominate over the line strengthening 
in such an MHD snapshot because of the systematic evacuation 
of magnetic elements whose interiors are expected to be consid¬ 
erably cooler than their hot surrounding walls at equal geometric 
height. 

4.1. MHD simulations versus idealized models 

The situation in a more realistic plage atmosphere is far more 
complex than the simple, heavily idealized models used in the 
past. This leads to the surprising result that horizontal RT not 
just weakens lines in the magnetic features, in some locations 
(mainly in the centers of magnetic elements), it actually strength¬ 
ens the lines. 

In idealized models, such as the thin-tube model, the wall 
may be defined in a simple way as the region where the transi¬ 
tion from the outer, hot, dense and field-free to the inner, cool, 
thin and strongly magnetized atmosphere takes place. In Paper I 
we defined the wall to be the region where the density decreases 
from 90% to 10% of the density outside the fiux-tube (i.e. of 
the non-magnetic quiet-Sun atmosphere at the same geometrical 
height), which also coincides with the magnetic field strength 
increase at the FT boundary. In the thin-tube approximation the 
temperature also decreases smoothly across the wall from very 
high values (10000 K and more) in the outer atmosphere towards 
low values (5000 K and less) in the inner parts of the FT/FS at 
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equal geometrical height. Therefore, the inner parts of the mag¬ 
netic element see hot walls. 

In the MHD snapshot the walls, i.e. the high-density re¬ 
gions from which the horizontal radiation reaching the core of 
the magnetic element originates, are generally cool, due to the 
presence of strong, cool down-flows around the magnetic ele¬ 
ments. The drop in density coincides with the increase of the 
magnetic held strength, but the temperature drop lies further 
away because the down-flows extend clearly beyond the mag¬ 
netic held. Therefore, the magnetic element is radiatively almost 
completely insulated from the outer hot atmosphere correspond¬ 
ing to granular upflows. Only little or no UV-irradiation from 
the hot granular gas reaches the core of sufficiently intense and 
large magnetic elements. This leads to a situation with strong 
UV irradiation from below, but a sub-average one from the side 
in the cores of magnetic elements. Consequently, the 1-D NLTE 
fines are weaker than their 3-D NLTE counterparts. That the UV 
irradiation from below is strong may also be seen from the EW 
ratios of the fines calculated in LTE and 1-D NLTE (see Eig.O. 
The 1-D NLTE lines are much weaker than those calculated in 
LTE as the latter ignore over-ionization due to the UV irradia¬ 
tion. Only near the boundary of the elements is the irradiation 
from the hot granules from the sides still effective and leads to 
the expected fine weakening. 

Because of the high density in the down-flows the continuum 
is formed at greater heights, i.e., lower temperatures in these 
walls of the magnetic element. This may also be seen in the con¬ 
tinuum image of Eig.[2l where the granules and the central part 
of ET/ES are brighter than the intervening lanes. 

In weaker magnetic elements, i.e. mainly in ES along inter¬ 
granular lanes, the line weakening is still the major effect, sim¬ 
ilar to the idealized models. Although there are cool and dense 
down-flows on two sides of the ES, the UV irradiation from the 
hot granule into the element still prevails. Towards the center, 
the under-population of Ee i levels due to horizontal UV irradi¬ 
ation is typically strongly suppressed, but a clear strengthening 
of the fines due to horizontal transfer is found only in a minority 
of cases. 

4.2. MHD versus HD model 

The plage model used in this work gives similar spatial averages 
for the strengths of the investigated lines as those we obtained 
in Paper II for the HD model with B = 0. Both, line weakening 
and strengthening, are present in each snapshot and, when aver¬ 
aged spatially over the whole snapshot, the differences between 
fines calculated in 1-D NLTE and those calculated in 3-D NLTE 
are small. On average, for every cool location in the atmosphere 
with a hot environment we also find a hot location with cool en¬ 
vironment. Thus, the global results for the atmospheres with a 
400 G average field and with no field are very similar, although 
they do not have the same origin when examined in detail. 

Eurthermore, in both models large differences may occur lo¬ 
cally, i.e. if spatially resolved spectra are considered. Local dif¬ 
ferences of several ten percent are found, those in the MHD 
snapshot being slightly stronger since the variation of atmo¬ 
spheric parameters, especially the density, is larger in the MHD 
model. The main differences between the two simulations occur 
in the intergranular lanes and their intersections. 

The magnetic field in the MHD model impedes the intergran¬ 
ular down-flows, which are therefore partially displaced towards 
the granules. The magnetic elements themselves are located in 
the middle of these cool and dense down-flows that separate the 
hot granular atmosphere from the magnetic element. This dimin¬ 


ishes the influence of the hot granular atmosphere on the interior 
of the magnetic element. In the core of a strong-held element 
it can even overcompensate for the systematic evacuation of the 
element. In the HD model, the transition from the granular to 
the intergranular atmosphere is smoother and resembles more 
the walls found in the classical thin-ET model. Therefore, in the 
HD model, the weakening of the lines computed in 3-D NLTE 
spans the whole intergranular lanes, though in a milder form, as 
the density difference between granules and intergranular lanes 
is not so large owing to the absence of a magnetic held. 

4.3. Are the FT/FS walls resolved? 

In Paper I we set the wall thickness of the mainly investigated 
ET model to be 20 km. Although, the situation in the MHD ET 
investigated here is more complex, we can still take the drop 
from 90 to 10 percent in density at the approximate geometrical 
height of Tc = 1 in the close granular environment as a measure 
for the wall thickness. Eor a ET with a horizontal extension of 
500 km and for a ES of similar width, the density drop typ¬ 
ically takes place over two to four voxels, i.e. over a distance 
of approximately 50 to 100 km (see Eigs.[8]and [TT]). Whether 
the relatively crude horizontal resolution of 21 km in the origi¬ 
nal MHD simulation manages to resolve the walls of magnetic 
elements is unclear. As the thickness of and the temperature dis¬ 
tribution in the wall is directly responsible for the strength of the 
effects of horizontal RT (see the results of the tests carried out 
in Paper I), it could well be that a better resolved snapshot may 
lead to a different quantitative result. We plan to test the influ¬ 
ence of resolution on the results of 3-D RT in MHD simulations 
in an upcoming investigation. 

5. Conclusions 

In this work, we continued our investigation of the influence 
of horizontal radiative transfer (RT) on diagnostically important 
Eei lines by calculating the line formation in a snapshot of a 
3-D radiation magneto-hydrody namic (MHD) simu lation per¬ 
formed with the MURAM code (IVogler et al.ll2005l) . The cho¬ 
sen snapshot represents a strong plage region (the simulation 
started with an initial vertical homogeneous held of 400 Gauss). 
By comparing profiles computed in LTE, 1-D NLTE and true 
3-D NLTE, we are able to isolate and quantify the effects 
(line weakening by UV irradiation from hot environments) of 
NLTE and horizon t al RT as described for the first time by 
IStenholm & Stenflol (Il977l) within a simple and highly idealized 
flux tube (ET) geometry. These r esults were confirmed within 
2-D (iBruls & von der Lfih3 l200lh and 3-D (Paper I) idealized 
models. We know from Paper II, that in pure hydrodynamic sim¬ 
ulations, the opposite effect, i.e. a strengthening of the 3-D NLTE 
fines takes place almost equally often, mainly inside granules. 
However, for a plage region one would naively expect a clear 
dominance of the line weakening effect due to the strong UV 
irradiation from hot walls into the evacuated magnetic elements. 

The surprising result of our computations is that the expected 
line weakening in the 3-D NLTE lines only takes place in the 
outer parts of a ET. At the center of the ET, the 3-D NLTE fines 
are either unaffected by the horizontal radiation inflow or are 
even strengthened. The reason for this unexpected behavior lies 
in the fact, that the ET in the MHD model are surrounded by 
dense and relatively cool down-flowing gas. This gas forms an 
isolating buffer between the magnetic element and the hot at¬ 
mosphere of the neighboring granules. The UV-irradiation from 
these hot parts of the photosphere, which dominates the behavior 
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in thin-tube models, is therefore strongly suppressed and has a 
noteworthy effect near the boundaries of the FT only. Even there, 
it is reduced due to the absorption by the intervening cool and 
dense down-flowing gas. The down-flows are usually not hotter 
than the inner part of the FT at equal geometrical height. The 
continuum in these dense down-flows is formed in higher and 
therefore cooler layers than in the granules or the innermost part 
of the FT. This can be clearly seen from the continuum image in 
Fig.El As the (UV-) continuum irradiation from these surround¬ 
ing down-flows is often lower than from directly below, the 3-D 
NLTE line is stronger than its 1-D NLTE counterpart which does 
not feel these ’’colder” contributions, from the side. 

In the case of flux sheets (FS) along the intergranular lanes, 
the situation is mostly similar. The magnetic field strength in 
the thinner FS is, however, often smaller 1 kG in a weak 
FS) and, therefore, the effects weaker. Line weakening due to 
horizontal irradiation takes place at the edges of the FS and in 
the surrounding dense lanes. The opposite effect in the center of 
the FS is weaker than for the stronger FT or only present in the 
form of a suppressed line weakening. 

In paper I we showed that the thickness of the down-flows 
and the wall, i.e. the region where the density drops and the mag¬ 
netic field increases when moving into the magnetic element, 
plays a central role for the qualitative result. The horizontal reso¬ 
lution in the employed MHD atmosphere is 21 km per voxel, i.e. 
much lower than that in the simplified models used in Paper I. In 
a future investigation, we plan to consider also simulation snap¬ 
shots with a finer grid. 

The significantly different behavior between the idealized 
thin-tube model analyzed in Paper I and the MHD simula¬ 
tions considered here introd uces a caveat for the uncritical 
use of the thin-tube model (iDefouwl Il976h . which has been 
widely and very successfully used. Besi des be i ng in valuable 
for studies of the thermal structure of (ISpruitl 119761) . MHD 
wave propagat ion in (e.,g. [Roberts & Webb Il979h. or, insta¬ 
bilities of (e .g. I Webb & Roberts! 1 19781: [Roberts & Wel^ll978l: 
ISpruitl[T97% magnetic elements, the thin-tube has turned out to 
describe a number of observations of polarized spec tra rather 
well (e.g., IZaver et ^ Il989l: iBruls & Solankil 1 19951) . Our re¬ 
sults indicate that the thin-tube approximation cannot capture 
the complexity of a magnetic element, so that at least in some 
contexts it should not be employed. However, the thin-tube ap¬ 
proximation provides a very good approximation of the mag- 
netic structure of magnetic e lements in 3-D MHD simulations 
(lYelles Chaouche et all1200^. at least if the second order thin- 
tube m odel is used fe.g. iPneuman et alii 19861: IFerriz-Mas'^ 
Il989l and later publications of these authors). This implies that 
it is mainly the structure of the external atmosphere which is re¬ 
sponsible for the difference between the present results and those 
of Paper I. 
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Fig. 1. Atmospheric parameters as indicated above each panel at heights z = 0 km (right panels) corresponding to < >= 1, 

and, z = 240 km, the approximate average line core formation region of the investigated lines (left panels). Positive velocities 
(v^ > 0 km) correspond to up-flows (blue). 
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Fig. 3. Spatial distribution of (left panels) and (right panels) for the 525.02 nm (upper panels) and the 630.15 nm 

line (lower panels). Positive indicates locations with weaker 3-D NLTE than the corresponding 1-D lines. 
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Fig. 4. Distribution of the (dashed lines) and (solid lines) values for the 525.02 nm (left panels) and the 630.15 nm 

line (right panels). The three vertically arranged panels contain the histograms for different spatial selections. Top panels: full 
atmosphere; Middle panels: Areas with strong up-flows (thin lines) and down-flows (thick lines), the spatial points were selected if 
the vertical velocity at = 1 was larger than 3 km/s in either direction. Bottom panels: The same as in the middle panels, but with 
the thick line now representing points at which the magnetic held strength \Bz\ is larger than 1 kG, and the thin line held strengths 
below 100 G. The number of selected spatial points in the middle and bottom panels amounts to approximately 10% to 20% of the 
whole atmosphere for each selection. 
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Fig. 5. Sample profiles at selected positions along the cut through the MHD FT. The x-positions from where the profiles originate 
are indicated above each panel. Line styles: thick solid: 3-D NLTE; thin solid: 1-D NLTE; dashed: LTE. Note the different y-axis 
ranges used in the panels. 
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Fig. 10. Temperature profile (left panel), vertical magnetic field (middle panel), and vertical velocity (right panel) in a vertical 
cut through a sample FS. The position of the cut is marked in the continuum image of Fig. [2] (vertical cut near the right upper edge). 
Additionally, the formation heights of the same spectral line cores as in Fig. [ 7 ] are given: Left and middle panel: Continuum (white 
line), 525.02 nm line core (yellow line), 630.15 nm line core (brown line). In the right panel the corresponding colors are chosen 
differently: white (continuum), blue (525.02 nm line core), and, black (630.15 nm line core). 
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